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ABSTRACT 

Three  "bef ore-and-af ter " studies  were  initiated  by  the 
State  of  Montana  Water  Quality  Bureau  in  1983  to  determine 
whether  the  upgrading  of  wastewater  treatment  facilities  at  Hot 
Springs,  Lewistown,  and  Hamilton  resulted  in  measurable 
improvements  in  the  water  quality  and  stream  biota  of  the 
receiving  streams  (Hot  Springs  Creek,  Big  Spring  Creek,  and  the 
Bitterroot  River) . Commonly,  the  effectiveness  of  improved 
wastewater  treatment  has  been  measured  only  by  the  improvement  in 
the  quality  of  the  discharge. 

Water  quality  data  and  benthic  macroinvertebrates  were 
collected  in  spring,  summer,  and  autumn  from  the  receiving  waters 
upstream  and  downstream  of  the  discharges  before  upgrading  began 
and  twelve  to  eighteen  months  after  upgrading  was  completed. 
Discharges  were  also  sampled  at  the  same  time.  The  limited 
amount  of  water  quality  data  precluded  a detailed  analysis  of 
changes  in  the  variables  measured.  Analysis  was  also  limited  by 
the  availability  of  flow  data  or  by  seasonal  and  year-to-year 
variability  of  streamflows.  Macroinvertebrates  were  analyzed 
using  a bioassessment  technique  which  integrates  eight  community- 
level  parameters . 

Improvement  in  the  quality  of  the  discharges  was 
measured  readily.  Generally,  BOD,  TSS,  and  ammonia 
concentrations  were  reduced  significantly.  Nitrate-nitrite 
concentrations  increased  significantly  because  of  more  complete 
treatment  of  influent  organic  material  and  more  thorough 
nitrification  of  ammonia.  Phosphorous  concentrations  decreased 
at  most  50%.  Improvement  in  the  quality  of  receiving  stream 
waters  was  less  apparent.  In  Hot  Springs  Creek,  variable 
streamflows  and  pollutant  loadings  from  other  sources  masked 
effects  of  the  Hot  Springs  wastewater  treatment  lagoon  discharge. 
Decreases  in  pollutant  concentrations  in  Big  Spring  Creek  up  to 
five  miles  downstream  of  the  Lewistown  WWTP  discharge  were 
discernable  but  subtle.  In  the  Bitterroot  River,  improvement 
could  be  detected  in  the  mixing  zone  but  two  miles  downstream  of 
the  Hamilton  WWTP  discharge  nutrient  concentrations  appeared  to 
increase . 


Analysis  of  macroinvertebrate  samples  detected,  at 
best,  only  modest  improvement  attributable  to  the  upgrading  of 
sewage  treatment  facilities.  Hot  Springs  Creek  and  Big  Spring 
Creek  appear  to  be  impaired  by  multiple  environmental  stresses 
which  largely  overwhelmed  the  improvement  attributable  to  higher 
quality  discharges.  The  impact  of  the  Hamilton  WWTP  on  the 
Bitterroot  River  was  slight  both  before  and  after  upgrading. 
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1.  INTRODUCTION 

Since  1972  nearly  200  million  dollars  have  been  spent 
in  Montana  under  the  Construction  Grants  Program  to  upgrade 
wastewater  treatment  plants  so  that  they  conform  with  secondary 
treatment  standards  required  by  the  Federal  Clean  Water  Act 
(Bahls  1982).  The  effectiveness  of  treatment  facilities  and 
upgraded  wastewater  treatment  is  often  judged  solely  on  whether  a 
facility  meets  its  effluent  limitations  (EPA  1985).  Improved 
instream  water  quality  and  uses  are  two  other  measures  of  the 
effectiveness  of  upgraded  wastewater  treatment.  In  Montana, 
little  or  no  attempt  had  been  made  to  document  these 
improvements . 

In  1983,  the  State  of  Montana  Water  Quality  Bureau 
(MWQB)  initiated  three  "bef ore-and-af ter " studies  to  document 
improvements  in  water  quality  and  stream  biota  resulting  from  the 
upgrading  of  wastewater  treatment  facilities  from  primary  to 
secondary  treatment.  The  three  wastewater  treatment  facilities 
selected  for  these  studies  serve  the  towns  of  Hot  Springs, 
Lewistown,  and  Hamilton.  These  facilities  were  selected  because 
the  ratio  of  discharge  flow  to  receiving  water  flow  was  larger 
than  it  was  for  other  facilities  scheduled  for  upgrading  at  the 
time.  Water  and  benthic  macroinvertebrate  samples  were 
collected  from  the  receiving  streams  in  1983  before  upgrading 
began  and  in  1986  and  1987  after  construction  was  completed. 
This  report  presents  the  results  of  sampling  conducted  on  Hot 
Springs  Creek,  Big  Spring  Creek  and  the  Bitterroot  River  - the 
receiving  waters  for  the  discharges  from  the  Hot  Springs, 
Lewistown  and  Hamilton  wastewater  treatment  facilities, 
respectively . 
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2.  OBJECTIVE 

The  objective  of  the  "bef ore-and-af ter " studies  was  to 
determine  whether  the  upgrading  of  wastewater  treatment 
facilities  at  Hot  Springs,  Lewistown  and  Hamilton  from  primary  to 
secondary  treatment  resulted  in  a measureable  improvement  in 
ambient  water  quality  and  stream  biota. 
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3.  STUDY  DESIGN 

Each  study  required  the  collection  of  water  quality  and 
benthic  macroinvertebrate  samples  before  upgrading  began  and 
after  upgrading  was  completed.  Samples  were  collected  12  to  18 
months  after  upgrading  to  allow  for  fine  tuning  of  treatment 
plant  operations  and  for  complete  recovery  of  the  aquatic 
ecosystem.  Monitoring  locations  included  one  control  site 
upstream  of  the  discharge,  one  station  just  below  the  discharge, 
and  one  or  two  stations  further  downstream.  Wastewater 
discharges  were  also  sampled.  Sampling  was  conducted  in  early 
spring,  late  summer,  and  mid-autumn. 

Upgrading  from  primary  to  secondary  treatment 
significantly  reduces  concentrations  of  total  suspended  solids 
( TSS ) and  biochemical  oxygen  demand  (BOD)  in  the  effluent. 
Secondary  treatment  may  also  reduce  potentially  toxic 
concentrations  of  ammonia  through  nitrification.  There  is  often 
no  significant  reduction  in  concentrations  of  nutrients 
(phosphorous  and  inorganic  nitrogen)  that  promote  algae  growth 
(Bahls  1982).  Water  samples  were  analyzed  for  nutrients  (ammonia 
and  nitrate-nitrite,  orthophosphorous , and  total  phosphorous), 
temperature,  pH,  total  suspended  solids  (TSS),  turbidity,  and 
specific  conductance.  Total  inorganic  nitrogen  (TIN)  was 
calculated  as  the  sum  of  the  nitrate-nitrite  and  ammonia 
concentrations.  Streamflow  and  discharge  flow  were  measured 
where  possible.  Self-monitoring  data  for  BOD  and  TSS  which  were 
collected  by  treatment  plant  operators  were  also  evaluated. 

Several  biotic  communities  were  analyzed.  Benthic 
macroinvertebrates  were  collected  because  the  structure  and 
composition  of  this  community  may  respond  to  a change  in  the 
quality  of  a discharge.  They  may  also  provide  a "truer" 
indication  of  "average"  ambient  water  quality  than  does  a limited 
water  sampling  program.  The  growth  of  periphyton  on  artificial 
substrates  in  Big  Spring  Creek  was  estimated  from  chlorophyll  & 
accrual  rates  to  determine  whether  changes  in  nutrient 
concentrations  in  the  Lewistown  WWTP  effluent  affected  primary 
production. 
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4 . METHODS 

4 . 1 Field  Sampling 

Water  temperature  was  measured  with  a mercury  field 
thermometer  having  gradations  of  1°C  and  a resolution  of  0.1°C. 
Streamflow  was  measured  with  a Marsh-McBirney  Model  2100  Water 
Current  Meter.  Wastewater  discharges  were  either  gaged  or 
determined  from  in-plant  flow  totalizers  or  weirs  and  rating 
tables . 

Benthic  macroinvertebrates  were  collected  with  a 
modified  Hess  sampler  for  quantitative  analysis,  or  with  a D-net 
using  a traveling  kick  or  unit  effort  kick  method  for  qualitative 
analysis.  Each  study  differed  in  the  method  used  and  number  of 
replicates  collected. 

Artificial  substrates  for  measuring  primary 
productivity  (Big  Spring  Creek  only)  consisted  of  eight  standard 
microscope  slides  held  in  commercial  plastic  carriages 
( Periphytometer  II)  which  were  anchored  to  the  substrate. 
Exposure  time  was  approximately  three  weeks.  Slides  were  then 
stored  frozen  until  chlorophyll  was  extracted  and  analyzed. 

Water  samples  were  collected  in  new  polyethylene 
bottles  and  preserved  following  guidelines  of  the  U.S. 
Environmental  Protection  Agency  (U.S.  EPA  1988). 

4.2  Water  Quality  Analysis 

The  analysis  of  water  quality  data  was  limited  by  the 
small  number  of  samples  which  were  collected.  At  most,  three 
samples  were  collected  at  each  site  before  and  after  upgrading 
in  the  Hot  Springs  Creek  and  Bitterroot  River  studies.  Five  or 
six  samples  were  collected  at  each  site  in  each  phase  of  the  Big 
Spring  Creek  study.  The  paucity  of  data  precluded  more  than  a 
non-statistical  comparison  of  yearly  averages  and  ranges  of 
concentrations . 

"Average"  concentrations  were  practically  meaningless 
in  Hot  Springs  Creek  because  of  extreme  seasonal  and  annual 
variation  in  flow,  and  effects  of  other  pollutant  sources. 
Concentrations  of  most  variables  fluctuated  wildly,  with  little 
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or  no  relationship  to  flow.  An  attempt  was  made  to  characterize 
"average"  water  quality  by  calculating  a standardized  loading. 
The  standardized  loading  describes  the  weight  of  a pollutant  in 
pounds  per  day  per  cubic  foot  per  second  (cfs)  of  water.  It  was 
determined  by  calculating  the  total  loading  of  a pollutant  (in 
pounds  per  day)  for  all  sampling  episodes  in  one  year  and 
dividing  it  by  the  total  flow  for  that  year.  The  standardized 
loading  is,  therefore,  nothing  more  than  a loading,  in  pounds  per 
cfs,  calculated  using  a flow-weighted  pollutant  concentration. 
The  same  result  would  be  achieved  by  physically  compositing  one 
second  of  flow  from  each  sampling  episode,  measuring  the 
pollutant  concentration  in  mg/l,  and  calculating  the  resulting 
load  per  cfs. 

4.3  Macroinvertebrate  Analysis 

Macroinvertebrate  data  were  evaluated  using  an 
integrated  analysis  technique  adapted  from  the  Rapid 
Bioassessment  Protocols  under  development  for  the  U.S. 
Environmental  Protection  Agency  (Plafkin  et  al.  1987).  Eight 
community-level  parameters  were  used  to  provide  a composite, 
integrated  evaluation  of  biologically  significant  degradation. 
These  parameters  were: 

1.  total  macroinvertebrate  abundance 

2.  taxa  richness 

3.  Shannon  diversity 

4 . biotic  index 

5.  EPT  richness  (number  of  distinct  taxa  in  the 
orders  Ephemeroptera,  Plecoptera,  and  Trichoptera) 

6.  ratio  of  EPT  to  Chironomidae  (relative  abundance 
of  Ephemeroptera,  Plecoptera,  and  Trichoptera  to 
Chironomidae ) 

7.  Percent  Community  Similarity  (PCS)  index 

8a.  Percent  Relative  Abundance  (PRA)  of  the  dominant 
taxa  (Hot  Springs  Creek  study  only) 
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8b.  ratio  of  scrapers  to  filter  feeders  (relative 
abundance  of  functional  feeding  groups,  Bitterroot 
River  and  Big  Spring  Creek  studies  only) 

These  parameters  provide  a wide  range  of  sensitivity  to 
water  quality  degradation.  Each  parameter  provides  a slightly 
different  measure  of  the  benthic  community  with  the  parameters' 
responses  varying  according  to  the  type  and  severity  of 
pollution.  The  sensitivity  of  some  parameters  overlap;  this 
redundancy  minimizes  the  consequences  of  an  aberrant  value  for 
one  parameter.  Integrating  individual  parameters  into  a 
composite  evaluation  increases  the  reliability  of  conclusions 
while  enhancing  the  sensitivity  of  the  evaluation  to  a wider 
range  of  pollution  effects. 

For  each  sampling  period  and  station,  all  parameters 
but  one  were  assigned  scores  based  on  their  similarity  to  the 
values  at  the  control  station  (Table  1)  . The  percent  relative 
abundance  of  the  dominant  taxa  was  assigned  a score  based  on  the 
actual  percent  contribution,  not  percent  comparability  to  the 
reference  station.  The  sum  of  the  scores  assigned  to  the  eight 
parameters  provides  an  integrated  evaluation  of  biologically 
significant  degradation  at  each  site  (Table  2). 

In  studies  where  replicate  samples  were  collected 
(Bitterroot  River  and  Big  Spring  Creek),  statistical  analyses 
were  used  to  test  for  differences  in  the  community-level 
parameters  among  sites  and  seasons.  These  tests  included  a two- 
way  Analysis  of  Variance  (Sorkal  and  Rohlf  1981)  and  the  Newman- 
Keuls  test  (Zar  1974). 

Where  appropriate,  total  macroinvertebrate  abundance, 
ordinal  densities,  ordinal  relative  abundances,  the  relative 
abundance  of  numerically  dominant  taxa  and  presence/absence  data 
were  compared  among  stations  and  between  years . 
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Table  1.  Criteria  for  assigning  scores  to  aquatic  macroinvertebrate 
parameters  based  on  their  percent  comparability  to  values 
obtained  from  the  control  station  (adapted  from  Plafkin  et 
al  1987)  . 


Scoring 

Criteria^ 

Parameter 

6 

4 

2 

0 

Abundance^ 

50-150% 

35-50% 

20-35% 

<20% 

or 

or 

or 

150-200% 

200-250% 

<250% 

Taxa  richness 

>80% 

60-80% 

40-60% 

<40% 

Shannon  diversity 

>90% 

80-90% 

70-80% 

<70% 

Biotic  index 

>90% 

80-90% 

70-80% 

<70% 

EPT  richness 
Ratio  of  EPT  to 

>90% 

80-90% 

70-80% 

<70% 

Chironomidae 
Percent  Community 

>75% 

50-75% 

25-50% 

<25% 

Similarity  (PCS) 
index 

>80% 

60-80% 

40-60% 

<40% 

Percent  Relative 

Abundance  (PRA)  of 
dominant  taxa^ 
Ratio  of  scrapers 

<20% 

25-35% 

35-45% 

>45% 

to  filterers^ 

>50% 

35-50% 

20-35% 

<20% 

1 Values  falling  at  break  points  for  scoring  categories  may  be  assigned 
intermediate  scores. 

Total  macroinvertebrate  abundance  substantially  higher  (>150%)  or 
substantially  lower  (<35%)  than  values  at  the  control  station  is 
indicative  of  environmental  degradation. 

^ Hot  Springs  Creek  study  only. 

^ Big  Spring  Creek  and  Bitterroot  River  studies  only. 
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Table  2.  Criteria  for  the  assessment  of  biologically  significant 
environmental  degradation  (from  Plafkin  et  al  1987). 


Cumulative  Score 

Classification 

40-48 

Non- impaired 

26-38 

Slightly  impaired 

10-24 

Moderately  impaired 

<8 

Severely  impaired 
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5.  HOT  SPRINGS  CREEK  STUDY 
5 . 1 Overview 

The  town  of  Hot  Springs  is  located  approximately  ninety 
miles  northwest  of  Missoula.  It  is  situated  on  the  southeast 
edge  of  the  Cabinet  Mountains  in  the  drainage  of  the  Little 
Bitterroot  River. 

The  new  wastewater  treatment  facility  replaces  an 
activated  sludge  plant  built  in  1954.  This  plant  had  a design 
flow  of  approximately  100,000  gallons  per  day  (gpd)  and  served 
approximately  600  persons.  In  the  late  1970's,  rehabilitation  of 
the  town's  sewer  lines  was  attempted  in  order  to  reduce 
infiltration  and  inflow,  but  plant  inflow  actually  increased 
after  the  work  was  completed.  By  1980,  approximately  80%  of  the 
flow  to  the  plant  was  infiltration,  and  bypasses  of  raw  sewage  to 
Hot  Springs  Creek  became  common.  The  plant  was  also  unable  to 
meet  effluent  limits  because  of  hydraulic  overload  which  reduced 
detention  time  throughout  the  treatment  process.  BOD 
concentrations  exceeded  limits  approximately  10  percent  of  the 
time  and  TSS  exceeded  limits  34  percent  of  the  time  (Storch 
Corporation/Engineers  1982). 

Construction  of  the  new  wastewater  treatment  facility 
began  early  in  1985  and  was  substantially  complete  by  mid-1986. 
The  facility,  which  has  a design  flow  of  237,000  gpd,  consists 
of  two  aerated  primary  treatment  cells  and  one  aerated  secondary 
treatment  cell,  followed  by  chlorine  disinfection.  Detention 
time  is  approximately  25  days. 

Hot  Springs  Creek  is  the  receiving  stream  for  the 
town's  wastewater  discharge.  It  originates  in  the  Cabinet 
Mountains  west  of  town,  and  joins  the  Little  Bitterroot  River 
approximately  five  miles  east  of  town.  It  is  the  only  major 
creek  in  the  area.  The  stream  is  classified  A-Closed  (protective 
of  public  surface  water  supplies)  upstream  of  the  Hot  Springs 
water  supply  intake  and  C-3  from  the  Hot  Springs  water  supply 
intake  to  the  Little  Bitterroot  River,  including  the  reach 
downstream  of  and  just  upstream  from  the  Hot  Springs  lagoon 
discharge.  Waters  classified  C-3  are  suitable  for  recreation. 
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and  growth  and  propagation  of  warm  water  aquatic  life.  They  are 
marginal  for  potable,  agricultural,  and  industrial  uses. 

Within  the  immediate  area  of  the  wastewater  treatment 
lagoon,  water  quality  in  Hot  Springs  Creek  is  strongly  influenced 
by  soil  structure,  soil  composition,  and  land  use  practices.  The 
dominant  soil  types  are  silts  and  clays  deposited  as  glacial  lake 
sediments.  These  soils  are  highly  erosive  and  very  susceptible 
to  undercutting  and  slumping.  Damming  and  ponding  by  beavers, 
which  is  quite  common,  alter  natural  stream  flows  and  divert  much 
of  the  stream  onto  adjacent  agricultural  fields.  This  sheet  flow 
reenters  the  main  creek  channel  at  diffuse  points  downstream  of 
the  wastewater  treatment  facility.  These  influences  are  the 
primary  determinants  of  water  quality  in  Hot  Springs  Creek  and 
may  limit  the  discernment  of  impacts  of  the  wastewater  discharge. 

5.2  Sampling  Locations  and  Description 

Three  monitoring  stations  were  located  on  Hot  Springs 
Creek.  The  control  (HS01)  was  located  approximately  one-quarter 
mile  upstream  of  the  discharge  (HS02).  One  station  (HS03)  was 
located  approximately  one  hundred  yards  downstream  of  the 
discharge.  Another  station  (HS04)  was  located,  approximately  two 
miles  downstream  of  the  discharge  (Figure  1). 

Several  stations  had  to  be  relocated  after  upgrading 
because  the  effluent  outfall  was  moved  approximately  three- 
quarters  of  a mile  downstream  of  its  original  location.  In  1987, 
the  control  site  was  located  just  downstream  of  1983 's  station 
HS03.  The  location  of  HS04  was  unchanged. 

5 . 3 Results 

5.3.1  Water  Chemistry 

Sel  f -monitoring  data  for  BOD  and  TSS  (Table  3) 
demonstrate  an  improvement  in  the  quality  of  the  wastewater 
discharge.  The  average  monthly  BOD  concentration  decreased  50% 
The  average  monthly  TSS  concentration  decreased  approximately 
25%.  TSS  concentrations  in  effluent  samples  collected  by  MWQB 
before  and  after  upgrading  (Table  4)  were  similar  to 


11 


Figure  1.  Sampling  locations  for  the  Hot  Springs  Creek  before-and-af ter  study. 
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Table  4.  Water  quality  data  for  the  Hot  Springs  Wastewater  Treatment  Lagoons  discharge  and  Hot  Springs  Creek,  1983  and  1987. 
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TSS  concentration  was  calculated  using  less  than  recommended  volume  of  water  and  is  subject  to  significantly  increased  error. 


Sampling  Location  Flow  Temperature  Lab  Turbidity  TSS  NH^  NO2+NO3  TIN  Ortho-P  Total-P  Specific 

and  Date  cfs  °C  pH  ntu  mg/L  mg/L-N  mg/L  mg/L-N  mg/L  mg/L  Conductance 

umhos  at  25° 
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concentrations  reported  by  the  permittee.  High  concentrations  of 
TSS  in  the  effluent  of  lagoon  treatment  systems  are  often  caused 
by  heavy  blooms  of  algae  which  thrive  in  these  nutrient-rich 
environments . 

The  average  ammonia  concentration  decreased  57%  in  the 
effluent  after  upgrading,  from  7.7  mg/1  to  3.3  mg/1,  but  the 
ranges  of  concentrations  indicate  that  ammonia  may  still 
approach  levels  observed  before  upgrading.  The  average  nitrate- 
nitrite  concentration  increased  from  0.19  mg/1  to  0.98  mg/1.  The 
average  orthophosphorous  concentration  increased  from  1.80  mg/1 
to  2.38  mg/1,  and  the  average  total  phosphorus  concentration 
decreased  from  2.93  mg/1  to  2.55  mg/1.  However,  the  ranges  of 
phosphorous  concentrations  suggest  that  there  was  little  if  any 
change  after  upgrading. 

Water  quality  in  Hot  Springs  Creek  was  difficult  to 
interpret  because  it  was  so  variable.  Generally,  it  was  more 
variable  in  1987  than  in  1983,  although  for  most  water  quality 
variables  two  of  the  three  concentrations  measured  in  1987  were 
within  the  ranges  observed  in  1983  (Table  4).  Because  the  third 
concentration  was  much  higher,  average  concentrations  of  ammonia, 
nitrate-nitrite,  TIN,  orthophosphorous  and  total  phosphorous  were 
greater  in  1987  than  in  1983  at  all  stream  locations  (Table  5). 
In  most  instances,  this  third  concentration  was  measured  in 
November,  1987  when  streamflow  was  extremely  low.  Average 
turbidity  and  TSS  were  lower  in  1987  at  HS01  and  HS03,  but  higher 
at  HS04 . 

5.3.2  Benthic  Macroinvertebrates 

One  benthic  macroinvertebrate  sample  was  collected  at 
each  stream  location  in  the  spring,  summer,  and  autumn  before  and 
after  upgrading,  using  the  traveling  kick  method.  Taxa 
comprising  at  least  10%  of  the  macroinvertebrate  fauna  at 
stations  above  or  below  the  discharge  are  listed  in  Table  6. 
Community  composition  at  HS03  was  highly  variable  between  years 
(Table  6)  and  among  sampling  dates  (unpublished  data).  Dramatic 
changes  in  the  abundance  of  dominant  taxa  indicate  frequent 


Table  5.  Average  values  for  selected  water  quality  variables  for  the  Hot  Springs  Wastewater  Treatment  Lagoons  discharge  and 
Hot  Springs  Creek,  1983  and  1987. 
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Table  6.  The  number  of  individuals  collected  annually  (all  samples 

combined)  and  the  percent  relative  abundance  of  taxa  comprising 
at  least  10%  of  the  macroinvertebrate  fauna  at  stations  above  or 
below  the  Hot  Springs  Wastewater  Treatment  Lagoons  during  1983 
and  1987. 


Taxon  Station: 

19  8 

3 

1 9 

8 7 

Total 

Number 

HS01 

HS03 

HS04 

Total 

Number 

HS0 1 

HS03 

HS04 

Tubificidae* 

0 

2929 

13% 

28% 

31% 

Simulium  spp. 

1238 

2% 

<1% 

45% 

1640 

9% 

19% 

10% 

Micropsectra  sp. 

274 

5% 

<1% 

5% 

2337 

13% 

34% 

3% 

Cheumatopsyche  sp. 

1156 

34% 

<1% 

11% 

1354 

30% 

2% 

9% 

Hyallela  azteca 

560 

14% 

2% 

6% 

1886 

14% 

5% 

33% 

Gammarus  sp. 

2154 

10% 

63% 

6% 

102 

<1% 

<1% 

3% 

Di.ame.ga  sp. 

806 

3% 

9% 

18% 

171 

3% 

<1% 

2% 

Prostoia  fresametga 

340 

10% 

<1% 

0% 

207 

5% 

1% 

0% 

* Tubificidae  were  not  identified  or  enumerated  in  1983  samples. 
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environmental  perturbations . 

Most  community-level  parameters  exhibited  relatively 
low  mean  values  characteristic  of  a moderately  impaired 
environment  (Table  7).  Some  parameters  clearly  improved  from 
1983  to  1987  (taxa  richness,  Shannon  diversity,  PRA  of  dominant 
taxa,  PCS  index).  Cumulative  scores  (Tables  8 and  9)  indicate 
improvement  at  downstream  stations  each  season  in  1987.  HS03 
was  moderately  impaired  two  of  three  seasons  both  years,  severely 
impaired  in  November  1983  and  slightly  impaired  in  April  1987. 
The  severe  impairment  in  November  1983  may  have  resulted  from  a 
spill  into  a sewer  of  1800  gallons  of  heating  oil,  which  was 
recovered  from  the  creek  between  HS02  and  HS03.  HS04  was 
moderately  impaired  two  seasons  in  1983  (slightly  impaired  the 
third  season)  and  slightly  impaired  two  seasons  in  1987 
(moderately  impaired  the  third).  In  both  years,  conditions  at 
HS04  appeared  to  improve  from  spring  to  autumn  (McGuire  1988a) . 

5.4  Discussion 

The  improvement  in  the  quality  of  the  Hot  Springs 
wastewater  discharge  is  evident  in  the  reduced  concentrations  of 
BOD  and  TSS . Phosphorous  concentrations  in  1987  were  comparable 
to  those  measured  in  1983.  Ammonia  appears  to  have  decreased  in 
the  effluent,  although  the  concentration  measured  in  November 
1987  is  almost  as  high  as  concentrations  measured  before 
upgrading.  Higher  nitrate-nitrite  concentrations  were  probably 
due  to  nitrification  of  ammonia.  It  is  likely  that  violations  of 
BOD  and  TSS  limits,  as  well  as  bypasses  of  raw  sewage  to  Hot 
Springs  Creek,  have  been  eliminated.  Pollutant  loadings  to  Hot 
Springs  Creek  have  also  decreased  dramatically  (Table  10). 

The  variability  of  pollutant  concentrations  and 
streamflow  makes  interpretation  of  instream  water  quality  data 
difficult.  Average  concentrations  for  most  variables  tended  to 
be  higher  in  1987  than  in  1983,  often  by  two  or  three  fold.  The 
very  poor  water  quality  and  abnormally  low  flows  which  occurred 
in  November  1987  skewed  average  concentrations  considerably. 
Standardized  water  quality  (pounds  per  day  per  cfs)  was  very 
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Table  7.  Mean  values  of  parameters  characterizing  macroinvertebrate 
assemblages  above  and  below  the  Hot  Springs  WWTP  during  1983  and 
1987  (three  samples  per  year  at  each  station) . 


PARAMETER 

1 

.983 

19  8 

7 

Station:  HS01 

HS03 

HS04 

HS0 1 

HS03 

HS04 

Abundance 

840 

920 

900 

1,060 

1770 

1,150 

Taxa  richness 

30 

18 

20 

28 

25 

27 

Shannon  diversity 

2.7 

1.9* 

2.6 

2.9 

2.6 

3.0 

Biotic  index  (0-10) 

5.1 

4.9 

5.6 

5.9 

6.5 

7.0 

EPT  richness 

6 

3 

2 

3 

2 

2 

EPT/C  ratio 

18.0 

0.02 

3.18 

2.72 

0.10 

1.71 

%RA  dominant  taxon 

48 

60* 

38 

33 

37 

29 

PCS  index 

— 

24* 

40 

— 

44 

46 

* based  on  March  and 

August  data  only. 

Table  8.  Scores  assigned  to  diagnostic  macroinvertebrate 
stations  below  the  Hot  Springs  WWTP  during  1983 
date) . 

parameters 
(one  sample 

at 

per 

Station: 

HS03 

HS04 

PARAMETER 

Mar . 

Aug . 

Nov . 

Mar 

Aug. 

Nov . 

Abundance 

6 

4 

2 

4 

0 

4 

Taxa  richness 

4 

4 

0 

2 

4 

6 

Shannon  diversity 

4 

0 

0 

0 

6 

6 

Biotic  index 

3 

6 

0 

3 

6 

6 

EPT  richness 

2 

0 

0 

0 

0 

6 

EPT/C  ratio 

0 

0 

0 

2 

0 

0 

PCS  index 

0 

0 

0 

0 

0 

2 

%RA  dominant  taxon 

4 

0 

0 

0 

3 

6 

Total  score 

23 

14 

2 

11 

21 

36 

20 


Table  9 . Scores  assigned  to  diagnostic  macroinvertebrate  parameters  at 
stations  below  the  Hot  Springs  Wastewater  Treatment  Lagoons 
during  1987  (one  sample  per  season) . 


Parameter 

Station: 

HS03 

HS04 

Apr . 

Aug . 

Nov . 

Apr . 

Aug. 

Nov. 

Abundance 

6 

5 

2 

4 

6 

4 

Taxa  richness 

4 

6 

2 

4 

6 

6 

Shannon  diversity 

4 

6 

4 

4 

6 

6 

Biotic  index 

4 

4 

6 

2 

5 

4 

EPT  richness 

0 

0 

6 

0 

0 

6 

EPT/C  ratio 

2 

0 

0 

0 

4 

6 

PCS  index 

2 

0 

2 

0 

2 

2 

PRA  dominant  taxa 

5 

2 

0 

6 

4 

2 

Total  score 

27 

23 

22 

20 

35 

36 

Table  10.  Average  flow  in  cubic  feet  per  second  (cfs),  and  average  loading  (AL)  in  pounds  per  day  and 

standardized  loading  (SL)  in  pounds  per  day  per  cfs.  Hot  Springs  Wastewater  Treatment  Lagoons 
and  Hot  Springs  Creek. 
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similar  both  years  at  HS01.  At  HS03,  standardized  loadings  of 
ammonia  and  TIN  in  1987  were  10%  of  1983  loadings. 
Orthophosphorous  and  total  phosphorous  were  two-thirds  lower. 
Similar  improvement  occurred  at  HS04  for  all  parameters  except 
TSS,  which  more  than  doubled. 

A number  of  other  influences  affect  water  quality  in 
Hot  Springs  Creek.  Highly  erodible  silts  and  clays  derived  from 
glacial  lake  sediments  form  much  of  the  creekbed  and  probably 
result  in  high  natural  background  levels  of  TSS  and  turbidity. 
Also,  beaver  damming  is  very  common,  altering  natural  streamflows 
and  nutrient  assimilation  and  cycling  processes.  This  activity 
also  diverts  much  of  the  streamflow  onto  adjacent  agricultural 
land.  This  flow  reenters  the  creek  at  diffuse  points  downstream, 
carrying  significant  amounts  of  sediment  and  nutrients.  Average 
TSS  and  nitrate-nitrite  concentrations  were  higher  at  HS04  than 
at  HS03  both  years.  Standardized  loads  of  TSS  clearly  increased 
both  years  between  stations . Nutrient  loadings  were  more 
variable . 

The  integrated  analysis  of  macroinvertebrate  samples 
indicated  that,  overall,  the  environmental  integrity  of  Hot 
Springs  Creek  improved  only  moderately.  Biotic  index  values  at 
HS03  and  HS04  in  1987  (6.5  and  7.0,  respectively)  indicated 
fairly  poor  water  quality  and  significant  organic  pollution. 
Macroinvertebrate  assemblages  in  both  years  were  characteristic 
of  a moderately  impaired  environment  and  multiple  stresses.  The 
apparent  improvement  documented  between  1983  and  1987  at 
downstream  stations  is  attributable,  in  part,  to  degradation  of 
the  control  site.  At  HS01,  sharp  declines  in  taxa  richness. 
Shannon  diversity,  and  EPT  richness  from  spring  to  summer  was 
possibly  due  to  chlorine  toxicity  (a  total  chlorine  concentration 
of  0.06  mg/1  was  measured  in  the  field  during  the  summer  of  both 
years).  Significantly  lower  ammonia  and  BOD  concentrations  in 
the  effluent,  as  well  as  the  elimination  of  raw  sewage  bypasses 
to  the  creek,  may  have  also  contributed  to  the  apparent 
improvement  in  macroinvertebrate  assemblages. 

Variability  in  community  composition  at  HS03  both  years 
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indicated  frequent  environmental  perturbations.  Stresses  on  the 
biological  community  at  this  site  caused  by  chlorine,  ammonia, 
and  nutrient  loading  were  at  least  additive  and  probably 
multiplicative.  Improvement  observed  at  HS04  from  spring  to 
autumn  both  years  may  be  attributable  to  winter  nitrite  toxicity 
from  ammonia  nitrification.  Amelioration  of  toxic  conditions 
followed  by  biological  recovery  could  have  then  occurred  in 
summer  and  autumn  (McGuire  1988a) . 


/ 
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6 . BIG  SPRING  CREEK  STUDY 
6 . 1 Overview 

The  town  of  Lewistown  is  located  in  Fergus  County 
approximately  100  miles  southeast  of  Great  Falls.  It  is 
situated  in  rolling  foothill  country  in  a mountainous  area  of  the 
central  part  of  the  state. 

Lewistown' s original  wastewater  treatment  plant  was 
built  in  1939  . The  plant  was  expanded  in  1965  and  provided 
primary  sedimentation  and  disinfection  by  chlorination  before 
discharging  to  Big  Spring  Creek.  At  the  time  of  expansion,  the 
collection  system  outfall  line  was  replaced  to  exclude  large 
volumes  of  infiltration.  This  plant  provided  adequate  primary 
treatment  during  dry  weather,  but  peak  wet  weather  flows  exceeded 
the  4 million  gallons  per  day  capacity  of  the  system,  requiring  a 
partial  bypass  of  raw  wastewater  to  the  chlorine  contact  basin. 

Further  upgrading  of  the  facility  occurred  and  was 
completed  in  August,  1986.  Secondary  treatment  is  provided  by 
rotating  biological  contactors  and  disinfection  by  ultraviolet 
radiation  before  wastewater  is  discharged  to  Big  Spring  Creek. 
The  primary  clarifiers  and  anaerobic  digesters  from  the 
previously  constructed  facilities  were  rehabilitated  and  utilized 
in  the  1986  upgrade.  The  plant  currently  serves  an  estimated 
7100  persons  and  is  designed  for  an  average  daily  flow  of  2.83 
million  gallons  and  a peak  hourly  flow  of  9.0  mgd.  Sewer  lines 
were  rehabilitated  in  1979,  and  infiltration  was  reduced  by 
approximately  0.7  0 mgd,  on  average,  and  by  0.95  mgd  at  peak 
flows.  However,  infiltration  and  inflow  still  contribute 
significantly  to  the  plant  flow  (HKM  Associates  1981). 

Big  Spring  Creek  is  the  receiving  stream  for  the 
Lewistown  WWTP  discharge.  The  creek's  source  is  a spring  located 
about  six  miles  southeast  of  Lewistown.  This  water  is  believed 
to  originate  as  subsurface  drainage  from  carbonate  rocks  in  the 
Big  Snowy  Mountains.  A fish  hatchery  is  located  near  the  spring. 
Big  Spring  Creek  flows  northwest  from  its  source,  through  the 
city  of  Lewistown,  to  its  confluence  with  the  Judith  River. 


25 


Big  Spring  Creek  was  described  as  a "hard,  calcium- 
bicarbonate"  type  water  (MDHES  1975),  reflecting  its  origins  in 
the  carbonate  terrains  of  the  Big  Snowy  Mountains.  The  creek 
experiences  less  variability  in  water  temperature,  discharge,  and 
concentrations  of  most  chemical  constituents  than  many  surface 
streams  because  of  the  relatively  constant  nature  of  groundwater 
recharge.  Big  Spring  Creek  is  classified  B-2  from  a point 
upstream  of  the  Lewistown  WWTP  to  the  Judith  River.  Waters 
classified  B-2  are  suitable  for  potable,  recreational, 
agricultural,  and  industrial  uses,  and  for  growth  and  marginal 
propagation  of  cold  water  aquatic  life.  The  stream  is  considered 
a "blue-ribbon"  trout  stream,  and  supports  populations  of  brown 
and  rainbow  trout.  A periphyton  analysis  of  the  stream  performed 
in  1974  reported  the  trophic  status  to  be  enriched,  probably  due 
to  organic  and  nutrient  loading  from  the  Lewistown  WWTP,  the  fish 
hatchery,  and  septic  tank  systems  (MDHES  1975). 

6.2  Sampling  Locations  and  Description 

Figure  2 shows  locations  sampled  upstream  and 
downstream  of  the  WWTP  discharge  (SC02).  The  station  upstream  of 
the  discharge  was  the  control  (SC01).  Station  SC03  was  located 
approximately  two  hundred  yards  downstream  of  the  wastewater 
discharge.  Station  SC04  was  located  approximately  two  miles 
downstream  of  the  discharge,  and  station  SC05  was  located 
approximately  five  miles  downstream. 

6 . 3 Results 

6.3.1  Water  Chemistry 

Self-monitoring  data  (Table  11)  indicate  that  the 
monthly  average  concentration  of  effluent  BOD  decreased  77%,  from 
38  mg/1  to  13  mg/1,  after  upgrading.  TSS  decreased  45%,  from  20 
mg/1  to  11  mg/1.  Monitoring  conducted  by  MWQB  (Tables  12  and  13) 
also  documented  a decrease  in  TSS  and,  in  addition,  turbidity 
(73%) . 

The  greatest  improvement  occurred  in  ammonia.  The 
average  concentration  decreased  88%  from  5.0  mg/1  in  1983  to  0.6 
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Figure  2.  Sampling  locations  for  the  Big  Spring  Creek  before-and  after  study. 
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Table  13.  Average  values  for  selected  water  quality  variables  for  the  Lewistown  WWTP  discharge  and  Big  Spring  Creek, 
1983  and  1987. 
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mg/1  in  1987.  The  average  nitrate-nitrite  concentration 
increased  almost  six-fold,  a result  of  ammonia  nitrification. 
Orthophosphorous  increased  slightly,  while  total  phosphorous 
decreased  slightly. 

Water  quality  in  Big  Spring  Creek  at  SC01  was  very 
similar  both  years.  Concentrations  of  most  variables  changed 
less  than  50%  from  1983  averages.  However,  ammonia  increased 
100%  and  orthophosphorous  decreased  50%  in  1987. 

Water  quality  downstream  of  the  WWTP  was  strongly 
influenced  by  the  quality  of  the  discharge.  The  clearest 
improvement  occurred  at  SC03.  Changes  in  TSS  and  nutrient 
concentrations  coincided  with  changes  in  the  discharge.  The 
ammonia  concentration  decreased  the  most  (85%);  nitrate-nitrite 
increased  the  most  (24%).  Water  quality  at  SC04  also  improved, 
although  the  change  in  most  pollutant  concentrations  was  smaller 
than  at  SC03.  AT  SC05,  all  pollutants  exhibited  a decrease  in 
average  concentration  after  upgrading,  including  nitrate-nitrite 
and  orthophosphorous,  which  had  increased  at  SC03  and  SC04  (Table 
13)  and  in  the  discharge. 

6.3.2  Benthic  Macroinvertebrates 

Three  modified  Hess  samples  and  one  traveling  kick 
sample  were  collected  at  each  stream  site  each  season.  Hess 
samples  which  were  collected  in  1983  were  sorted  and  enumerated 
to  the  ordinal  level  only.  All  other  samples  were  sorted  and 
enumerated  to  genus  or  species. 

Data  from  the  kick  samples  were  used  in  the  integrated 
analysis  of  community-level  parameters  for  comparing  before  and 
after  conditions.  These  same  parameters  were  calculated  for  the 
1987  Hess  samples  and  were  analyzed  by  a two-way  Analysis  of 
Variance  for  significant  differences  by  location  or  season.  A 
Newman-Keul ' s test  was  used  to  compare  each  pair  of  factors 
within  that  variable. 

Based  on  the  1983  data,  slight  degradation  of  the 
aquatic  environment  was  discernible  at  each  station  below  the 
WWTP  (Table  14).  With  a single  exception,  cumulative  scores  were 
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Table  14.  Mean  values  of  parameters  characterizing  macroinvertebrate 
assemblages  in  Big  Spring  Creek  above  and  below  the  Lewistown 
WWTP  during  1983  and  1987  (three  kick  samples  per  station) . 


STATION 


PARAMETER 

SC01 

SC03 

SC04 

SC05 

Abundance 

1983 

2,027 

2,879 

2,657 

2,137 

1987 

2,713 

2,562 

4,360 

2,665 

Taxa  richness 

1983 

35 

36 

32 

27 

1987 

29 

31 

27 

27 

Shannon  diversity 

1983 

3.37 

2.96 

2.96 

2.76 

1987 

3.52 

3.45 

3.06 

3.51 

Biotic  index  (0-10) 

1983 

4.02 

5.14 

4.75 

4.50 

1987 

5.22 

5.36 

4.91 

4.90 

EPT  richness 

1983 

16 

14 

13 

12 

1987 

9 

10 

10 

10 

EPT-Chironomidae  ratio 

1983 

10.9 

4.3 

13.7 

10.4 

1987 

1.61 

0.94 

2.68 

3.26 

Scraper-f ilterer  ratio 

1983 

0.25 

0.41 

0.17 

0.31 

1987 

0.63 

0.52 

0.69 

0.73 

PCS  index 

1983 

— — — 

55 

47 

50 

1987 

— 

46 

33 

43 

Percent  community  similarity 

(X100) 

using  station 

SC01  as 

the  basis  for 

comparison . 
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in  the  slightly  impaired  category  (Table  15).  The  October  score 
for  station  SC03  was  slightly  lower  (25),  falling  between  the 
slightly  and  moderately  impaired  ranges.  Mean  scores  for 
stations  SC03,  SC04  and  SC05  were  31,  35  and  35  respectively. 
Overall,  EPT  richness  was  the  most  sensitive  parameter.  A steady 
decline  in  the  number  of  mayfly,  stonefly  and  caddisfly  species 
collected  from  SC01  to  SC05  (Table  14)  lowered  scores  for  each 
downstream  station.  The  percent  community  similarity  and  biotic 
indices  also  indicated  slight  impairment  at  each  downstream  site. 
The  EPT  to  Chironomidae  ratio  was  depressed  at  SC03. 

The  1987  data  indicated  improved  conditions  downstream 
of  the  WWTP  when  compared  to  the  upstream  control  site  (Table 
16).  Mean  scores  for  the  1987  data  were  40,  36,  43  at  stations 
SC03,  SC04  and  SC05,  respectively.  For  each  date  tested  during 
1987,  station  SC04  was  rated  as  slightly  impaired  while  stations 
SC03  and  SC05  were  in  the  non-impaired  category.  Overall,  scores 
for  six  of  the  eight  parameters  improved  between  1983  and  1987. 
EPT  richness  and  the  biotic  index  accounted  for  43%  and  29%, 
respectively,  of  the  net  increase.  Higher  values  assigned  to  the 
scraper  to  filterer  ratio  (functional  feeding  groups),  Shannon 
diversity,  taxa  richness,  and  EPT  to  Chironomidae  ratio  also 
contributed  to  the  higher  cumulative  scores . 

However,  there  is  evidence  that  degradation  of  the 
control  site  contributed  to  the  higher  scores  during  1987.  Taxa 
richness,  biotic  index,  EPT  richness  and  EPT  to  Chironomidae 
ratio  indicated  worsened  conditions  at  SC01  during  1987  compared 
to  1983  (Table  14).  Since  these  parameters  accounted  for  85%  of 
the  net  improvement  at  downstream  sites,  the  validity  of  the 
assessment  for  1987  is  questionable. 

To  test  the  hypothesis  that  environmental  conditions 
had  deteriorated  between  1983  and  1987,  the  1983  data  were  used 
as  controls  for  assigning  scores  to  each  kick  sample  collected 
during  1987  (Table  17).  Because  laboratory  methods  differed 
between  years  (e.g.  subsampling  in  1987;  the  inclusion  of 
tubificids  and,  apparently,  more  Chironomidae  in  the  analyses  for 
1987)  some  parameters  were  biased.  However,  methods  were 
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Table  15.  Scores  assigned  to  diagnostic  macroinvertebrate  parameters  at 


stations  below 
per  season) . 

the 

Lewistown 

WWTP 

during 

1983 

( one 

kick 

sample 

Station: 

Parameter 

SC03 

SC04 

SC05 

Apr 

. Aug . 

Oct 

. Apr . Aug . 

Oct . 

Apr . 

Aug . 

Oct . 

Abundance 

6 

0 

6 

6 

4 

6 

6 

6 

6 

Taxa  richness 

6 

6 

5 

6 

6 

6 

6 

4 

4 

Shannon  diversity 

4 

6 

4 

4 

2 

6 

4 

0 

6 

Biotic  index 

4 

2 

2 

4 

3 

4 

3 

4 

6 

EPT  richness 

4 

6 

0 

2 

4 

1 

4 

2 

0 

EPT/C  ratio 

2 

2 

2 

6 

6 

6 

4 

6 

6 

Scraper/f ilterer  ratio 

4 

6 

4 

6 

6 

2 

6 

6 

6 

PCS  index 

6 

4 

2 

4 

2 

2 

4 

2 

3 

Total  score 

36 

32 

25 

38 

33 

33 

37 

30 

37 

Mean  score 

31 

35 

35 

Table  16.  Scores  assigned 
stations  below 
per  season) . 

to 

the 

diagnostic  macroinvertebrate  parameters  at 
Lewistown  WWTP  during  1987  (one  kick  sample 

Station: 

Parameter 

SC03 

SC04 

SC05 

Mar.  Aug. 

Oct . 

Mar , 

. Aug . 

Oct . 

Mar 

. Aug . 

Oct . 

Abundance 

6 

6 

4 

4 

6 

4 

6 

4 

6 

Taxa  richness 

6 

6 

6 

6 

6 

6 

6 

6 

6 

Shannon  diversity 

6 

4 

6 

4 

0 

6 

4 

4 

6 

Biotic  index 

6 

6 

6 

4 

6 

6 

6 

6 

6 

EPT  richness 

2 

6 

6 

6 

6 

6 

6 

6 

6 

EPT/C  ratio 

4 

2 

4 

6 

6 

4 

6 

6 

6 

Scraper/f ilterer  ratio 

6 

6 

6 

4 

2 

6 

6 

6 

6 

PCS  index 

4 

4 

2 

4 

0 

0 

4 

2 

2 

Total  score 

40 

40 

40 

38 

32 

38 

44 

40 

44 

Mean  score 

40 

36 

43 

36 


Table  17.  Bioassessment  scores  for  stations  on  Big  Spring  Creek  above 

and  below  the  Lewistown  WWTP  during  1987  using  1983  data 
as  controls  (one  kick  sample  at  each  station  per  season) . 


Station: 

SC01 

SC03 

SC04 

SC05 

Mar 

. Aug . 

Oct . 

Mar . 

Aug . 

Oct . 

Mar . 

Aug . 

Oct . 

Mar . 

Aug . 

Oct . 

Abundance 

4 

6 

6 

6 

2 

4 

2 

6 

4 

4 

2 

4 

Taxa  richness 
Shannon 

6 

6 

4 

6 

4 

6 

4 

6 

6 

6 

6 

6 

diversity 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

Biotic  index 

1 

4 

3 

2 

6 

6 

4 

6 

6 

4 

4 

6 

EPT  richness 

0 

0 

0 

0 

4 

4 

2 

4 

6 

0 

4 

6 

EPT/C  ratio 

0 

4 

6 

0 

2 

6 

0 

6 

0 

0 

6 

2 

Scraper/ f ilterer 

6 

6 

6 

6 

6 

6 

4 

6 

6 

6 

6 

6 

ratio 
PSC  index 

3 

2 

2 

2 

4 

0 

4 

2 

2 

4 

2 

6 

Total  score 

26 

34 

33 

28 

34 

38 

26 

42 

36 

30 

36 

42 

Mean  score 

31 

33 

35 

36 

* 
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consistent  for  each  year  and  these  biases  should  have  been  the 
same  at  each  station.  Therefore,  the  mean  cumulative  score  for 
each  station  should  reflect  the  degree  of  degradation  between 
years,  with  lower  scores  indicating  relatively  more  impact  than 
higher  scores.  The  mean  score  at  SC01  (31)  was  lower  than  at 
progressive  downstream  stations  by  two,  four  and  five  points, 
respectively . 

If  this  analysis  is  valid  and  if  impacts  were  additive, 
the  influence  of  control  site  degradation  (reduced  standards  for 
comparison)  on  scores  at  each  station  between  1983  and  1987  can 
be  quantified.  Immediately  below  the  WWTP  effluent  (SC03),  the 
mean  score  increased  from  31  during  1983  to  40  during  1987  , so 
approximately  20%  (2/9ths)  of  this  apparent  improvement  resulted 
from  degradation  at  the  control  site.  Thus,  nearly  80%  of  the 
increased  score  was  attributable  to  improved  water  quality. 
Further  downstream,  improved  environmental  conditions  were  either 
not  detected  (SC04)  or  were  primarily  attributable  (approximately 
60%  at  SC05)  to  reduced  standards  of  comparison. 

Statistical  analysis  of  Hess  samples  collected  in  1987 
failed  to  provide  a clear  pattern  of  differences  among  stations 
(Table  18).  No  significant  differences  (p  = 0.05)  were  detected 
for  mean  Shannon  diversity  or  EPT  richness  among  stations.  The 
scraper  to  filterer  and  EPT  to  Chironomidae  ratios  were 
significantly  higher  at  SC05  than  at  other  stations.  These 
differences  were  due  to  fewer  midges  and  more  Helicopsyche  sp.  (a 
caddisfly  scraper)  at  the  downstream  station  and  probably 
indicate  more  stable  substrate  as  well  as  reduced  organic 
pollution.  At  station  SC04,  total  macroinvertebrate  abundance 
was  significantly  higher  while  the  ratio  of  insects  to  tubificids 
was  significantly  lower  than  at  other  stations.  These  results 
indicate  that  organic  enrichment  from  the  WWTP  effluent  had  its 
greatest  impact  on  benthic  invertebrates  at  this  site. 
Stoneflies  and  mayflies  were  also  more  abundant  at  this  site, 
suggesting  the  primary  effect  of  the  effluent  at  SC04  was  to 
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Table  18.  Synopsis  of  statistical  analyses  (ANOVA  and  Newman-Keuls  tests) 
of  selected  parameters  characterizing  macroinvertebrate 
assemblages  at  stations  above  and  below  the  Lewistown  WWTP 
during  1987  (nine  modified  Hess  samples  per  station)1. 


Parameter 

Significant  differences  among  locations 

Abundance 
Taxa  richness 
Shannon  diversity 
Biotic  index 
EPT  richness 
EPT-Chironomidae  ratio 
Scraper-f ilterer  ratio 
Insect-Tubif icidae  ratio 

SC01,  SC03 , SC05  < SC04 
SC05  < SC03 

no  significant  differences 
SC01  < SC03 , SC04 , SC05 
no  significant  differences 
SC01,  SCO 3 , SC04  < SC05 
SC01,  SC03 , SC04  < SC05 
SC01,  SC03 , SC05  > SC04 

The  March,  1987  Hess  sample  #1  from  HS04  was  destroyed  during  transport. 
The  pooled  (when  appropriate)  or  mean  parameter  values  of  the  remaining 
replicates  were  substituted  for  this  missing  value  for  statistical 
analyses . 
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enhance  the  food  supply.  The  preponderance  of  tubificids  and 
chironomids  at  SC04  was  partially  attributable  to  the  finer 
substrate  in  this  stream  reach. 

The  biotic  index  provided  a reliable  measure  of  organic 
enrichment  attributable  to  the  WWTP.  The  mean  biotic  index  value 
at  SC01  (4.04)  was  significantly  lower  than  at  the  downstream 
stations  and,  based  on  Hilsenhoff’s  (1987)  scale  for  evaluation, 
was  indicative  of  very  good  water  quality  with  possible  slight 
organic  pollution.  Mean  values  at  SC03,  SC04  and  SC05  were 
5.11,  5.96  and  5.06,  respectively.  Values  between  4.51  and  5.50 
indicate  good  water  quality  with  some  organic  pollution  while 
values  between  5.51  and  6.50  typify  fair  water  quality  and  fairly 
significant  organic  pollution  (Hilsenhoff  1987). 

Quantitative  comparisons  of  species  abundance  between 
years  were  impeded  by  differences  in  laboratory  processing  and 
taxonomy.  Organisms  in  the  1983  Hess  samples  were  not  identified 
below  order,  leaving  the  unreplicated  kick  samples  as  the  only 
source  of  data  for  discrete  taxa  from  both  years.  Discrepencies 
between  these  data  sets  further  complicated  efforts  to  interpret 
the  data.  Tubificidae  ( Oligochaeta ) , which  during  1987  comprised 
one,  twenty-two,  fifteen  and  five  percent  of  the  fauna  in  kick 
samples  at  stations  SC01,  SC03,  SC04  and  SC05,  respectively,  were 
not  included  in  the  1983  data.  Therefore,  tubificids  were 
excluded  when  calculating  the  relative  abundances  in  Table  19. 
For  most  numerically  dominant  taxa,  site  specific  changes  in 
mean  relative  abundance  between  1983  and  1987  were  within  the 
range  of  variability  normally  expected.  The  seven-fold  increase 
in  the  relative  abundance  of  Physa  sp.  at  station  SC01  was  an 
exception.  However,  kick  samples  may  not  have  provided  a 
reliable  measure  of  the  snails'  relative  abundance.  While  Physa 
sp.  densities  increased  at  SC01  between  1983  and  1987,  they  were 
relatively  less  abundant  in  Hess  samples  (three  percent)  than  in 
kick  samples  (23%). 

Areawide  differences  in  the  abundance  of  some  species 
were  more  distinct  (Table  19)  and  were  not  contradicted  by  the 
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Table  19.  Percent  relative  abundance  of  numerically  dominant  macro- 
invertebrate taxa  in  traveling-kick  samples  from  Big  Spring  Creek 
above  and  below  the  Lewis  town  WWTP  during  1983  and  1987  (three 
samples  per  station  each  year  combined) . 


1 9 

8 3 

1 9 

8 7 

Taxon 

SC01 

SC03 

SC04 

SC05 

SC01 

SC03 

SC04 

SC05 

DIPTERA 

QrthQc.lacii.ua  spp. 

9% 

11% 

2% 

2% 

3% 

3% 

5% 

2% 

Polypedilum  sp. 

4% 

11% 

9% 

5% 

1% 

5% 

1% 

1% 

Antocha  sp. 

6% 

2% 

<1% 

<1% 

7% 

1% 

0% 

<1% 

EPHEMEROPTERA 

Baetis  tricaudatus 

6% 

17% 

21% 

33% 

12% 

9% 

20% 

29% 

Ephemerella 
inf requens / inermis 

12% 

17% 

21% 

7% 

10% 

7% 

24% 

13% 

COLEOPTERA 

Optioserviis 

quadrimaculatus 

6% 

1% 

1% 

5% 

5% 

2% 

1% 

5% 

PLECOPTERA 

.Lsoper.la 

quinquepunctata 

10% 

4% 

5% 

4% 

<1% 

<1% 

4% 

1% 

TRICHOPTERA 

Hydropsyche 

occidentalis? 

18% 

5% 

10% 

23% 

7% 

2% 

1% 

13% 

Brachycentrus 

americanus 

14% 

4% 

3% 

1% 

4% 

2% 

2% 

<1% 

GASTROPODA 
Physa  sp. 

3% 

11% 

8% 

2% 

23% 

8% 

0% 

2% 

ISOPODA 

Asellus  sp. 

0% 

<1% 

5% 

3% 

1% 

12% 

9% 

• 4% 

TURBELLARIA 

1% 

2% 

2% 

4% 

3% 

17% 

10% 

6% 
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1987  Hess  data.  The  relative  abundance  of  one  stonefly  ( Isoperla 
quequinpunctata ) and  two  caddisflies  (Hydropsyche  occidentalis 
and  Brachycentrus  americanus ) declined  at  all  four  stations  in 
1987  compared  to  1983  while  the  converse  was  true  for  flatworms 
(Turbellaria)  and  an  isopod  ( Asellus  sp.).  These  changes  may 
indicate  habitat  degradation  throughout  the  study  area.  During 
both  years,  several  abundant  species  appeared  to  be  affected  by 
the  WWTP  effluent.  The  abundance  of  the  cranefly  Antocha  sp., 
and  the  riffle  beetle  Optioservus  quadrimaculatus  f were  reduced 
downstream  of  the  WWTP,  while  the  midge  Polypedilum  sp.  increased 
immediately  downstream  of  the  WWTP  effluent. 

Comparing  species  lists  for  the  two  years  provided 
further  evidence  of  environmental  deterioration  at  SC01  between 
1983  and  1987.  Twelve  species  of  stoneflies,  mayflies  and 
caddisflies  ( Alloperla  sp.,  Claassinia  sp.,  Cultus  sp.,  Isoperla 
(possibly  pints) , Baetjs  quilleri , Nixe  simplicoides , Heptagenia 
soltflri,  Pars  1 eptQp.h  1 eb  ia  (probably  heteronia)  , Symphitopsyche 
cockerelli , &.  slossonae , Brachycentrus  occidentalis  f and 
Glossosoma  sp.)  present  at  SC01  in  1983  were  not  found  in  1987 
(either  kick  or  Hess  samples).  All  twelve  species  were  rare 
during  1983  with  only  one  or  two  individuals  of  most  species 
being  collected.  Consequently,  some  species  may  have  been  missed 
due  to  subsampling.  However,  Cultus  sp.,  Paraleptophlebia  sp. 
and  Glossosoma  sp.  were  present  at  downstream  stations  during 
1987  and  only  one  species  of  caddisfly  (Hydroptila  sp.)  was 
collected  at  SC01  in  1987  that  was  not  found  there  in  1983.  The 
loss  of  so  many  sensitive  species  must  be  considered  evidence  of 
increased  environmental  stress  at  the  control  site. 

Data  from  Hess  samples  also  suggested  diminished 
environmental  quality  in  1987  relative  to  1983.  Mean 
macroinvertebrate  densities  declined  by  41%,  65%  and  60%  at 
stations  SC01,  SC03,  and  SC05,  respectively,  between  years  (Table 
20).  At  SC04,  where  the  mean  macroinvertebrate  density  increased 
three  percent  between  1983  and  1987,  tubificids  replaced 
caddisflies  as  the  predominant  benthic  invertebrate.  Caddisfly 
densities  and  relative  abundances  were  greatly  reduced  at  all 
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Table  20.  The  mean  density  (number  per  sample)  and  percent  relative 
abundance  (in  parentheses)  of  macroinvertebrate  orders  in  Big 
Spring  Creek  above  and  below  the  Lewistown  WWTP  during  1983  and 
1987  (nine  modified  Hess  samples  per  station  each  year) . 


STATION 


TAXON 


SC01  SC03  SC04  SC05 


Plecoptera 


1983 

52 

(3) 

61 

(3) 

24 

(2) 

14 

(1) 

1987 

8 

(1) 

2 

(<1) 

59 

(4) 

21 

(4) 

Trichoptera 

1983 

902 

(57) 

522 

(24) 

704 

(46) 

883 

(58) 

1987 

353 

(38) 

142 

(19) 

159 

(10) 

196 

(32) 

Ephemeroptera 

1983 

199 

(13) 

523 

(24) 

214 

(14) 

350 

(23) 

1987 

126 

(13) 

80 

(10) 

174 

(11) 

139 

(23) 

Coleoptera 

1983 

40 

(3) 

13 

(1) 

13 

(1) 

53 

(4) 

1987 

100 

(11) 

36 

(5) 

10 

(1) 

36 

(6) 

Diptera 

1983 

358 

(23) 

831 

(38) 

253 

(17) 

144 

(10) 

1987 

295 

(32) 

236 

(31) 

359 

(23) 

71 

(12) 

Isopoda 

1983 

0 

1 

(<1) 

149 

(10) 

16 

(1) 

1987 

1 

(<1) 

48 

(6) 

143 

(9) 

96 

(16) 

Gastropoda 

1983 

5 

(<1) 

181 

(8) 

34 

(2) 

10 

(1) 

1987 

27 

(3) 

29 

(4) 

6 

(<1) 

1 

(<1) 

Other 

1983 

13 

(1) 

36 

(2) 

50 

(3) 

43 

(3) 

1987 

24 

(3) 

187 

(25) 

661 

(42) 

48 

(8) 

Total 

1983 

1987 

1570 

934 

2166 

760 

1516 

1571 

1514 

608 
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stations  during  1987.  Mayfly  densities  were  similarly  reduced 
during  1987  and  stoneflies  were  less  abundant  at  SC01  and  SC03  in 
1987.  Apparent  changes  in  relative  abundance  should  be  evaluated 
cautiously  because  Tubificidae  were  included  in  the  1987  data 
(McGuire  1988b) . 

6.3.3  Chlorophyll  £ Production 

The  rate  of  growth  of  algal  biomass  on  artificial 
substrates  was  estimated  by  measuring  the  rate  of  production  of 
chlorophyll  Chlorophyll  a generally  comprises  one  to  two 
percent  of  algal  biomass.  Estimates  of  algal  biomass  production 
were  derived  by  multiplying  the  chlorophyll  a accrual  rate  by  a 
constant  of  67,  which  represents  1.5%  of  algal  biomass  as 
chlorophyll  a (Table  21) . 

The  greatest,  and  only  significant  (p=.05)  difference 
in  chlorophyll  a accrual  between  1983  and  1987  occurred  at  SC03, 
where  the  average  accrual  rate  and  estimated  biomass  production 
were  almost  four  times  greater  in  1987.  Average  values  for  both 
variables  were  slightly  higher  at  SC01,  and  lower  at  SC04  and 
SC05,  in  1987,  but  not  significantly  so. 

6.4  Discussion 

Concentrations  of  BOD,  TSS,  and  ammonia  in  the  effluent 
decreased  significantly  after  upgrading.  Nitrate-nitrite 
concentrations  increased  almost  six-fold,  a result  of  ammonia 
nitrification.  Concentrations  of  orthophosphorous  and  total 
phosphorous  were  only  slightly  changed. 

Water  quality  in  Big  Spring  Creek  reflected  the 
improvement  in  the  quality  of  the  discharge  (Table  13).  At  SC03, 
water  quality  was  strongly  affected  by  the  improvement  in  the 
quality  of  the  discharge.  At  SC04  and  SC05,  the  WWTP  discharge 
was  still  a major  determinant  of  water  quality,  in  that  nutrient 
concentrations  had  not  decreased  to  background  levels.  However, 
nutrient  concentrations  at  these  stations  in  1987  were  usually 
slightly  less  than  or  unchanged  from  concentrations  in  1983. 


Table  21.  Chlorophyll  a accrual  and  estimated  biomass  production,  Big  Spring  Creek,  1983  and  1987. 
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Improvement  in  instream  water  quality  could  not  be  so 
easily  and  incontrovertibly  documented  by  macroinvertebrate 
populations.  An  assessment  of  the  affects  of  improved  water 
quality  was  complicated  by  altered  conditions  at  the  control 
site . 

Prior  to  upgrading,  the  WWTP  effluent  had  a slight  but 
easily  identifiable  impact  on  the  biota  of  the  receiving  stream. 
Some  sensitive  species  experienced  localized  reductions  in 
abundance  near  the  effluent  outfall  (e.g.  Opt ioservus 
quadrimaculatus  and  Hydropsyche  occidentalis ) while  others  were 
less  abundant  at  all  downstream  stations  (e.g.  Antocha  sp.  and 
Brachycentrus  americanus ) . Tolerant  species  were  more  abundant 
immediately  below  the  discharge  (e.g.  Polypedilum  sp.  ) or 
throughout  the  downstream  study  area  (e.g.  Baetis  tricaudatus  and 
Ephemerella  inermis /inf requens ) . In  addition,  many  community- 
level  indices  (total  macroinvertebrate  density.  Shannon 
diversity,  biotic  index  and  EPT  to  Chironomidae  ratio)  exhibited 
classic  responses  to  mild  organic  pollution  at  the  station 
immediately  downstream  of  the  WWTP's  discharge.  Further 
downstream,  the  slight  impairment  detected  was,  at  least 
partially,  a consequence  of  poor  habitat. 

During  1987,  the  abundance  and  taxa  richness  of  aquatic 
insects  was  reduced  at  all  stations  in  Big  Spring  Creek  relative 
to  1983,  including  SC01.  Insect  densities  were  reduced  by  40%  to 
60%  throughout  the  study  area.  The  most  dramatic  impact  was 
upstream  of  the  WWTP  effluent,  where  twelve  species  of 
stoneflies,  mayflies  or  caddisflies  present  during  1983  were  not 
found  in  1987.  Habitat  requirements  of  the  species  most  affected 
indicated  sediment  deposition  as  the  probable  cause  of  the 
significant  environmental  deterioration  reported  for  1987.  The 
severity  of  this  impact  may  have  been  augmented  by  the  lack  of 
periodic  flushing  provided  by  high  stream  flows  and  the 
precipitation  of  large  quantities  of  calcite  (calcium  carbonate). 
Sediment  deposition  and  calcite  precipitation  appeared  to  be  more 
severe  in  1987  than  in  1983.  Calcite  precipitation  appeared  to 
be  greatest  at  SC01,  and  less  noticeable  or  not  present  below  the 
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outfall.  Sediment  deposition  and  substrate  shifting  and  scouring 
was  evident  throughout  the  study  reach. 

The  integrated  bioassessment,  although  compromised  by 
lower  standards  of  comparison  resulting  from  degradation  at  the 
control  site,  detected  no  impairment  at  two  of  the  three 
downstream  stations . The  improved  sewage  treatment  did  appear  to 
reduce  environmental  stress  in  downstream  reaches,  particularly 
near  the  effluent  outfall  (SC03).  The  biotic  index  indicated 
some  organic  pollution  downstream  of  the  WWTP  during  1987.  The 
effects  of  organic  enrichment  were  most  pronounced  at  station 
SC04  approximately  two  miles  downstream.  For  most  numerically 
dominant  taxa,  abundance  and  distribution  patterns  were  similar 
between  1983  and  1987,  although  differences  among  stations  were 
generally  diminished.  A slight  impact  attributable  to  organic 
enrichment  from  the  WWTP  was  present  during  1987;  however,  this 
impact  was  largely  overwhelmed  by  deterioration  of  the  physical 
environment  in  Big  Spring  Creek  (McGuire  1988b) . 

In  1987,  chlorophyll  a accrual  just  downstream  of  the 
outfall  was  almost  four  times  higher  than  before  upgrading. 
Average  accrual  rates  at  SC04  and  SC05  were  slightly  lower, 
though  not  significantly  so.  The  increase  in  chlorophyll  & 
accrual  and  estimated  biomass  production  at  SC03  may  have  been 
due  to  higher  concentrations  of  bio-available  nutrients  in  the 
discharge,  especially  during  the  testing  period  (March  and 
April).  Ammonia  nitrif ication  increased  nitrate-nitrite 
concentrations  in  the  effluent  almost  ninefold  after  upgrading. 
Orthophosphorous  increased  slightly.  This  increase  was  barely 
detectable  at  SC03  possibly  because  of  rapid  nutrient 
assimilation  by  algae  and  macrophyte  growth  in  this  reach.  The 
relatively  constant  nutrient  concentrations  at  SC03,  SC04,  and 
SC05  may  indicate  that  the  assimilative  capacity  of  the  stream 
had  been  exceeded  or  reached. 
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7.  BITTERROOT  RIVER  STUDY 

7 . 1 Overview 

Hamilton's  WWTP , which  serves  a town  of  population  2700 
located  in  the  Bitterroot  River  valley  about  50  miles  south  of 
Missoula,  was  built  in  the  1950's.  Excessive  infiltration  and 
hydraulic  overloading  necessitated  reconstruction  of  the  WWTP  to 
meet  a fourfold  increase  in  influent  volume  that  occurs  during 
the  summer  months . 

The  new  WWTP  is  an  oxidation  type  activated  sludge 
facility.  Two  concentric  oxidation  ditches  and  three  secondary 
clarifiers  provide  the  treatment  capacity  to  handle  peak  summer 
flows.  As  a result  of  summer  irrigation,  infiltration  increases 
influent  flows  from  0.7  mgd  in  the  winter  to  2.5  mgd  in  the 
summer.  After  final  clarification,  the  effluent  passes  through  a 
chlorine  contact  basin  before  discharging  to  the  Bitterroot 
River.  Dechlorination  of  the  effluent  by  sulphur  dioxide  is 
practiced  in  the  summer. 

The  Bitterroot  River  is  a relatively  large  stream  that 
generally  provides  a dilution  ratio  of  at  least  100:1  for 
Hamilton's  wastewater  discharge.  During  spring  runoff,  dilution 
may  reach  1000:  1.  Due  to  the  linear  nature  of  the  stream  reach 
just  downstream  of  the  WWTP  outfall,  the  effluent  plume  hugs  the 
right  bank  for  upwards  of  one  mile  before  the  channel 
configuration  promotes  significant  lateral  mixing. 

The  river's  headwaters  lay  in  the  Bitterroot  Range 
south  of  Hamilton.  The  river  flows  northward  to  Missoula  where 
it  joins  with  the  Clark  Fork  River.  The  Bitterroot  River  is 
classified  B-l.  It's  waters  are  suitable  for  potable, 
recreational,  agricultural,  and  industrial  uses,  and  for  growth 
and  propagation  of  cold  water  aquatic  life.  The  river  and  valley 
are  known  for  their  scenic  and  recreational  opportunities  with 
hunting,  fishing,  camping,  and  boating  being  very  popular. 

7.2  Sampling  Locations  and  Description 

Three  monitoring  stations  were  located  on  the 
Bitterroot  River.  The  control  (BR01)  was  located  approximately 
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one-quarter  mile  upstream  of  the  discharge  (BR02).  The  station 
approximately  one-half  mile  downstream  of  the  discharge  (BROS) 
was  located  within  the  effluent  plume.  Conditions  at  tnis 
location  reflect  the  worst  possible  before  upgrading,  with  the 
most  potential  for  improvement  after  upgrading.  The  farthest 
downstream  station  (BR04)  was  approximately  two  miles  below  the 
discharge  (Figure  3). 

7 . 3 Results 

7.3.1  Water  Chemistry 

Water  quality  data  (Tables  22  and  23)  indicate  that 
turbidity,  TSS,  and  ammonia  were  reduced  over  90%  after 
upgrading.  Orthophosphorous  and  total  phosphorous  were 
approximately  half  their  concentrations  before  upgrading. 
Nitrate-nitrite  concentrations,  which  were  barely  detectable 
prior  to  upgrading,  increased  almost  two  hundred  times.  TIN 
increased  almost  fifty  percent,  possibly  because  of  more  complete 
treatment  of  organic  material  previously  discharged  as  TSS. 

Loadings  of  TSS,  ammonia,  orthophosphorous,  and  total 
phosphorous  to  the  Bitterroot  River  (Table  24)  decreased 
significantly.  Nitrate-nitrite  loading  increased  almost  seventy- 
five  fold.  Increases  in  TIN  and  decreases  in  phosphorous  and  TSS 
loadings  indicate  more  complete  nitrification  and  removal  of 
organic  material. 

Limited  water  quality  data,  and  the  lack  of  flow  data 
for  the  Bitterroot  River,  precluded  a detailed  analysis  of 
instream  water  quality  conditions.  Generally,  water  quality  at 
BR03,  which  was  located  in  the  plume  of  the  effluent,  reflected 
the  quality  of  the  discharge.  Ammonia,  orthophosphorus,  and 
total  phosphorous  decreased  commensurately  with  their  effluent 
concentrations . Nitrate-nitrite  concentrations  increased 
sixfold.  TSS  and  turbidity  were  influenced  more  by  water  quality 
upstream  of  the  discharge  than  by  the  improvement  in  the  quality 
of  the  discharge. 

Nutrient  concentrations,  excepting  ortho-phosphorous, 
were  slightly  higher  in  1986  than  1983  at  BR04 . Turbidity  and 
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Figure  3.  Sampling  locations  for  the  Bitterroot  River  before  and-after  suudy. 
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Table  24. 


TSS  and  nutrient  loading  (pounds  per  day)  to  the  Bitterroot 
River  from  the  Hamilton  WWTP,  1983  and  1986. 1 


DATE 

TSS 

nh3-n 

N02+N03-N 

TIN 

0-P04 

T-P 

March  1983 

274 

58.9 

0.1 

59.0 

39.4 

45.3 

August  1983 

437 

45.0 

2.9 

47.9 

46.6 

67.5 

November  1983 

224 

42.1 

0.2 

42.3 

28.9 

49.1 

Average 

312 

48.7 

1.0 

49.7 

38.3 

54.0 

April  1986 

7 

0.3 

69.7 

70.0 

19.6 

21 . 1 

August  1986 

— 

7.4 

63.0 

70.4 

22.9 

30.2 

November  1986 

17 

0.4 

107.4 

107.8 

31.9 

33.6 

Average 

12 

2.7 

80.0 

82.7 

24.8 

28.3 

Percent  Change 

-96 

-94 

+7520 

+ 66 

-35 

-85 

1 Loadings  were 

calculated  using 

instantaneous 

flows 

measured  at 

the 

time  of  sampling. 
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TSS,  as  at  BR03,  were  most  affected  by  ambient  water  quality 
upstream  of  the  discharge. 

7.3.2  Benthic  Macroinvertebrates 

In  1983,  one  composite  kick  sample  was  collected  per 
site  each  season.  In  1986,  three  unit-effort  kick  samples  were 
collected  per  site  each  season. 

During  1983,  the  mean  values  of  four  community-level 
parameters  indicated  differences  in  macroinvertebrate  assemblages 
above  and  below  the  WWTP  outfall  (Table  25).  Some  of  these 
indices  (biotic  index,  EPT  to  Chironomidae  ratio  and  PCS  index) 
returned  to  ambient  levels  further  downstream  (BR04).  Based  on 
Hilsenhoff's  (1987)  evaluation  of  water  quality  using  biotic 
index  values,  stations  BR01  and  BR04  had  excellent  water  quality 
with  no  apparent  organic  pollution  while  BR03  had  very  good  water 
quality  with  possible  slight  organic  pollution.  The  relative 
increase  in  the  abundance  of  filter  feeders  (expressed  as  the 
ratio  of  scrapers  to  filter  feeders)  suggested  that  increased 
nutrient  loading  changed  the  foodbase  downstream  from  the  WWTP. 

Population  data  for  numerically  dominant  taxa  also 
indicated  an  impact  attributable  to  the  WWTP  effluent  during  1983 
(Table  26).  The  relative  abundance  of  Lepidostoma  sp.,  the  most 
numerous  macroinvertebrate  in  the  study  area,  was  severely 
reduced  at  BR03  compared  with  the  other  stations.  Conversely, 
Polypedilum  sp.,  a species  more  tolerant  of  organic  enrichment, 
was  a minor  component  of  the  macroinvertebrate  assemblage  at  BR01 
and  BR04,  but  was  the  second  most  abundant  species  at  BR03  where 
it  comprised  13%  of  the  fauna.  The  virtual  elimination  of 
Lepidostoma  sp.  at  BR03  may  have  been  due  to  toxic  conditions 
(probably  episodes  of  high  chloramines  concentrations)  near  the 
effluent  outfall.  Lepidostoma  may,  therefore,  be  useful  for 
detecting  chlorine  plus  ammonia  impacts  in  streams  where  it  is 
abundant . 

Following  the  WWTP  upgrade,  differences  in  the  relative 
abundances  of  numerically  dominant  taxa  were  not  evident  (Table 
26)  although  values  of  community-level  parameters  varied  among 
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Table  25.  Mean  values  of  parameters  characterizing  macroinvertebrate 
assemblages  above  and  below  the  Hamilton  WWTP  during  1983  and 
1986  (three  samples  per  station  in  1983  and  nine  samples  per 
station  in  1986 ) . 


1983  1986 


PARAMETER  Station: 

BR01 

BR03 

BR04 

BR01 

BR03 

BR04 

Abundance 

3,300 

2,500 

3,600 

2,800 

2,900 

4,800 

Taxa  richness 

48 

47 

44 

58 

51 

49 

Shannon  diversity 

3.4 

3.5 

3.1 

3.7 

3.7 

3.4 

Biotic  index  (0-10) 

2.9 

4 . 1 

2.7 

3.1 

3.4 

3.6 

EPT  richness 

28 

27 

27 

30 

28 

23 

EPT/C  ratio 

21 

10 

21 

8 

4 

4 

Scraper/f ilterer  ratio 

2.1 

0.8 

0.6 

3.9 

1.2 

1.3 

PCS  index  * 

55 

65 

71 

69 

Percent  community  similarity 

(X100) 

using 

station  BR01 

as  the 

basis  for 

comparison. 


Table  26.  The  mean  percent  relative  abundance  of  the  numerically  dominant 

macroinvertebrate  taxa  above  and  below  the  Hamilton  WWTP  during 
1983  and  1986  (three  samples  per  station  in  1983  and  nine 
samples  per  station  in  1986)  and  for  all  samples  combined. 


1983  1986 


Taxon 

Overall 

BR01 

BR03 

BR04 

BR01 

BR03 

BR04 

Lepidostoma  sp. 
Ephemerella 

29% 

32% 

2% 

44% 

31% 

37% 

32% 

infrequens 

15% 

21% 

25% 

18% 

8% 

6% 

15% 

Optioservus  spp. 
Hydropsyche 

7% 

10% 

10% 

7% 

7% 

2% 

5% 

slossonae  ? 
Brachycentrus 

4% 

3% 

6% 

4% 

5% 

2% 

4% 

occidentalis 

4% 

4% 

3% 

12% 

1% 

1% 

1% 

Polypedilum  sp. 

3% 

<1% 

13% 

2% 

1% 

2% 

3% 

Drunella  flavilinea 

3% 

0% 

0% 

0% 

4% 

9% 

4% 

Hydroptila  sp. 

2% 

<1% 

0% 

1% 

7% 

5% 

2% 

Baetig  tricaudatus 

2% 

2% 

6% 

1% 

1% 

2% 

1% 
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stations  (Table  25).  The  two-way  analysis  of  variance  on 
replicate  samples  indicated  significant  differences  (p  = 0.05) 
among  stations  for  five  of  the  seven  parameters  while  four 
parameters  exhibited  significant  seasonal  variation.  The  Newman- 
Keuls  test  identified  significant  differences  among  stations  for 
only  three  parameters:  total  abundance,  biotic  index  and  EPT  to 
Chironomidae  ratio  (Table  27).  The  multiple  range  test  did  not 
identify  significant  differences  among  stations  for  mean  Shannon 
diversity  or  EPT  richness  values  although  the  ANOVA  found 
significant  differences  in  both  parameters  (p  = 0.05).  The 
inconsistent  statistical  results  were  probably  a consequence  of 
skewed  data  distributions. 

Using  the  1983  data  as  a baseline,  the  1983  and  1986 
data  for  the  control  station  were  compared  to  determine  their 
similarity  and  to  test  for  upstream  perturbations  that 
influenced  the  study  reach.  Scores  for  the  1986  data  ranged 
from  44  to  34  with  a mean  of  41.  The  spring  and  autumn  values, 
as  well  as  the  overall  mean,  were  in  the  non-impaired  category 
which  indicated  the  data  sets  were  generally  comparable. 
However,  the  cumulative  score  for  the  summer  sample  was  in  the 
slightly  impaired  range.  The  lower  score  appeared  attributable 
to  slight  ambiguities  in  field  and  laboratory  procedures  rather 
than  to  altered  environmental  conditions.  For  instance,  smaller 
and  therefore  more  midges  were  apparently  included  in  the  1986 
data.  While  anomalies  between  data  sets  precluded  more  detailed 
comparisons  between  the  1983  and  1986  data,  they  did  not  bias  the 
comparative  analysis. 

Based  on  the  1983  data,  a slight  impairment  of  the 
aquatic  environment  was  discernable  at  both  stations  below  the 
WWTP  (Table  28).  The  impact  was  most  evident  at  BR03  where,  on 
each  occasion,  the  cumulative  score  was  28,  near  the  bottom  of 
the  slightly  impaired  range.  Scores  were  more  variable  at  BR04, 
ranging  from  30  to  42.  The  November  sample  from  BR04  received  a 
non-impacted  rating  while  the  spring  and  summer  scores  were  in 
the  slightly  impaired  category. 


* 
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Table  27.  Synopsis  of  statistical  analyses  ( Newman-Keuls  test)  of 
selected  parameters  characterizing  macroinvertebrate 
assemblages  at  stations  above  and  below  the  Hamilton 
WWTP  during  1986. 


Parameter 

Location 

Abundance 

BR01  = BR03  < BR04 

Taxa  richness 

BR01  = BR03  = BR04 

Shannon  diversity 

BR01  = BR03  = BR04 

Biotic  index 

BR01  < BR04 

EPT  richness 

BR01  = BR03  = BR04 

EPT/C  ratio 

BR01  > BR03  = BR04 

Scraper/f ilterer  ratio 

BRO 1 = BR03  = BR04 

Table  28.  Scores  assigned  to  diagnostic  macroinvertebrate  parameters  at 
stations  below  the  Hamilton  WWTP  during  1983  (one  sample  per 
season) . 


Station:  BR03 


BR04 


PARAMETER 


Mar . Aug . Nov . Mar . Aug . Nov . 


Abundance 
Taxa  richness 
Shannon  diversity 
Biotic  index 
EPT  richness 
EPT/C  ratio 

Scraper/f ilterer  ratio 
PCS  index 


6 

6 

6 

0 

4 

2 

2 

2 


2 

6 

6 

2 

6 

2 

2 

2 


4 

4 

6 

0 

2 

4 

6 

2 


6 

6 

6 

4 

6 

6 

0 

4 


0 

6 

4 

6 

6 

6 

0 

2 


6 

6 

4 

6 

6 

4 

6 

4 


Total  score 
Mean 


28  28  28 
28 


38  30  42 

37 
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Following  the  WWTP  upgrade,  scores  for  BR03  were 
consistently  higher  (Table  29)  though  still  in  the  slightly 
impaired  category  on  two  of  the  three  dates  sampled  in  1986. 
With  the  mean  score  increasing  from  28  in  1983  to  37  in  1986, 
environmental  stresses  appear  to  have  been  reduced  near  the  WWTP 
outfall.  At  station  BR04,  approximately  2 miles  downstream,  the 
mean  cumulative  score  decreased  from  37  in  1983  to  31  in  1986  . 
This  decline,  while  not  significant  with  regard  to  impact 
severity,  was  in  accord  with  the  slight  increase  in  nutrient 
concentrations  measured  at  the  downstream  station  during  1986 
compared  to  1983.  Thus,  the  upgraded  treatment  at  the  Hamilton 
WWTP  appears  to  have  reduced  impacts  near  the  outfall  while 
shifting  the  impact  associated  with  nutrient  loading  further 
downstream  (McGuire  1988c). 

7.4  Discussion 

Water  quality  data  clearly  indicate  that  the  WWTP 
upgrading  improved  the  quality  of  the  discharge.  TSS  was  reduced 
over  90%.  Phosphorous  concentrations  and  loadings  to  the 
Bitterroot  River  also  decreased.  Inorganic  nitrogen  loading 
increased,  probably  because  of  more  complete  treatment  of 
influent  organic  nitrogen.  However,  the  ammonia  nitrogen 
fraction  decreased  to  less  than  0.5  ppm  because  of  more  complete 
nitrification . 

An  improvement  in  instream  water  quality  was  best 
documented  at  BR03.  As  this  station  was  located  in  the  effluent 
plume,  changes  in  nutrient  concentrations  were  consistent  with 
their  increase  or  decrease  in  the  WWTP  effluent.  Two  miles  below 
the  discharge  (BR04),  average  concentrations  of  all  nutrients, 
except  orthophosphorous , increased  in  1986.  This  increase  was 
especially  perplexing  because  ammonia,  TIN,  and  total 
phosphorous  had  decreased  at  BR03.  However,  nitrate-nitrite, 
orthophosphorous,  and  total  phosphorous  did  exhibit  a station-to- 
station  decrease  between  BR03  and  BR04.  No  satisfactory 
explanation  could  be  developed  to  explain  these  results.  One 
possibility  is  that  a different  benthic  community  became 
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Table  29.  Scores  assigned  to  diagnostic  macroinvertebrate  parameters  at 
stations  below  the  Hamilton  WWTP  during  1986  (three  samples  per 
season) . 


Station:  BR03 


BR04 


PARAMETER 


Apr . Aug . Nov . Apr . Aug . Nov . 


Abundance 
Taxa  richness 
Shannon  diversity 
Biotic  index 
EPT  richness 
EPT/C  ratio 

Scraper/f ilterer  ratio 
PCS  index 


4 

6 

6 

6 

6 

4 

4 

6 


6 

6 

6 

6 

6 

6 

0 

4 


6 

4 

6 

4 

2 

2 

6 

4 


0 

6 

6 

4 

4 

2 

2 

6 


4 

6 

6 

4 

4 

2 

2 

4 


6 

4 

6 

4 

0 

4 

4 

2 


Total  score 
Mean 


42  36  34 

37 


30  32  30 

31 
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established  in  response  to  changes  in  nitrogen  and  phosphorous 
concentrations,  which  resulted  in  different  nutrient 
assimilation  dynamics.  The  net  effect,  nonetheless,  was  a 
slightly  more  enriched  environment  at  BR04  after  upgrading  than 
before . 

Results  of  the  integrated  bioassessment  of 
macroinvertebrate  communities  suggest  that  the  upgraded  treatment 
appears  to  have  reduced  impacts  near  the  outfall  while  shifting 
the  impact  associated  with  nutrient  loading  further  downstream. 
Impacts  identified  at  BR03  in  1983  included  a cumulative  score 
near  the  bottom  of  the  slightly  impaired  range,  compared  to  a 
slightly  impaired  or  non-impaired  rating  at  BRQ4;  a biotic  index 
value  indicative  of  very  good  water  quality  with  possible  slight 
organic  pollution  compared  to  excellent  water  quality  with  no 
apparent  organic  pollution  at  BR04;  a severe  localized  reduction 
in  the  caddisfly  Lepidostoma  sp.  , possibly  from  chloramine 
toxicity;  and  a localized  abundance  of  the  midge  Polypedilum  sp., 
which  is  tolerant  to  organic  pollution.  A less  distinct  impact 
was  detected  at  station  BR04  in  the  spring  and  summer  of  1983, 
where  enhanced  periphyton  growth  affected  the  community 
structure . 

In  1986,  cumulative  scores  were  consistently  higher  at 
BR03  but  lower  at  BR04,  although  both  were  slightly  impaired. 
The  lack  of  apparent  differences  in  the  relative  abundances  of 
numerically  dominant  taxa  suggest  that  the  impact  at  BR04  was 
relatively  slight  compared  to  the  impact  at  BR03  in  1983. 
Increased  scraper/f ilterer  ratios  suggest  less  effect  from 
nutrient  and  organic  loading. 

In  1983,  impacts  below  the  WWTP  discharge  were  most 
likely  related  to  organic  enrichment  and  chloramine  or  ammonia 
toxicity.  In  1986,  impacts  associated  with  effluent  toxicity 
were  no  longer  detectable.  Impacts  associated  with  nutrients  and 
organic  enrichment,  while  slight  both  years,  were  more  evident  at 
BR04  in  19  86  . Overall,  the  Bitterroot  River  appears  to  be  in 
excellent  condition  with  macroinvertebrate  assemblages  indicative 
of  aquatic  environments  with  very  good  habitat  and  water  quality. 
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8.  CONCLUSIONS 

In  each  before-and-af ter  study,  an  improvement  in  the 
quality  of  the  WWTP  discharge  attributable  to  upgrading  could  be 
readily  documented.  In  all  three  cases,  TSS  loadings  were 
reduced.  Ammonia  concentrations,  often  high  enough  to  be 
potentially  toxic,  were  reduced  by  nitrification.  This  resulted 
in  much  higher  loadings  of  nitrate-nitrite  nitrogen  to  the 
receiving  streams. 

Documenting  improvement  in  the  quality  of  the  receiving 
streams  was  more  difficult.  In  Hot  Springs  Creek,  extremely 
variable  streamflows  and  significant  flow  alterations  and 
pollutant  loadings  from  other  sources  masked  the  effects  of  the 
lagoon  discharge.  In  Big  Spring  Creek,  decreases  in  pollutant 
concentrations  could  be  detected  at  three  stations,  and  up  to 
five  miles,  downstream  of  the  discharge.  Improvement  was  more 
subtle  as  downstream  distance  increased.  In  the  Bitterroot  River 
a significant  improvement  could  be  detected  in  the  mixing  zone, 
but  two  miles  downstream  nutrient  concentrations  may  have 
actually  increased  after  upgrading. 

Analysis  of  macroinvertebrate  samples  detected,  at 
best,  only  modest  improvement  attributable  to  the  upgrading  of 
sewage  treatment  facilities.  Hot  Springs  Creek,  before  and  after 
upgrading,  was  impaired  by  multiple  environmental  stresses, 
including  habitat  degradation  and  organic  pollution  related  to 
natural  biological  and  geological  processes.  Under  a more  normal 
flow  regime  results  of  the  upgrading  may  be  more  apparent. 
Improved  sewage  treatment  at  Lewistown  appeared  to  reduce 
environmental  stress  in  downstream  reaches;  however,  this 
improvement  was  largely  overwhelmed  by  deterioration  of  the 
physical  environment  throughout  the  study  reach.  Impacts  of  the 
Hamilton  WWTP  on  the  Bitterroot  River  appeared  to  be  slight  both 
before  and  after  upgrading.  As  a result  of  more  complete 
nitrification  after  upgrading,  impacts  associated  with  nutrient 
enrichment  were  shifted  downstream.  Overall,  the  river  was  and 
continues  to  remain  in  excellent  condition. 
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